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We consider the problem of discriminating two different quantum states in the setting of asymptotically many 
copies, and determine the optimal strategy that minimizes the total probability of error. This leads to the identifi- 
cation of the quantum Chernoff bound, thereby solving a long standing open problem. The bound reduces to the 
classical Chernoff bound when the quantum states under consideration commute. The quantum Chernoff bound 
is the natural symmetric distance measure between quantum states because of its clear operational meaning and 
because of the fact that it does not seem to share the undesirable features of other distance measures like the 
fidelity, the trace norm and the relative entropy. 
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INTRODUCTION 

One of the most basic tasks in information theory is the dis- 
crimination of two different probability distributions: given a 
source that outputs variables following one out of two possi- 
ble probability distributions, determine which one it is with 
the minimal possible error. In a seminal paper, Chernoff (J] 
solved this problem in the asymptotic regime and showed that 
the probability of error P e in discriminating two probability 
distributions decreases exponentially in the number of tests n 
that one can perform: P e ~ exp(— n£cs)- The optimal ex- 
ponent £c b arising in the asymptotic limit is called the Cher- 
noff bound [2]. One of the virtues of the Chernoff bound is 
that it yields a very natural distance measure between proba- 
bility distributions; it is essentially the unique distance mea- 
sure in the ubiquitous situation of independent and identically- 
distributed (i.i.d.) random variables. 

A quantum generalization of this result is highly desired. 
Indeed, the concept of randomness is much more elementary 
in the field of quantum mechanics than in classical physics. 
Given the large amount of experimental effort in the context 
of quantum information processing to prepare and measure 
quantum states, it is of fundamental importance to have a the- 
ory that allows to discriminate different quantum states. De- 
spite considerable effort, this quantum generalization of the 
Chernoff bound has until now remained unsolved. The prob- 
lem is to discriminate two sources that output many identical 
copies of one out of two different quantum states p and a, and 



the question is to identify the exponent arising asymptotically 
when performing the optimal test to discriminate them. This 
task is so fundamental that it was probably the first problem 
ever considered in the field of quantum information theory; it 
was solved in the one-copy case more than 30 years ago |yl|4|]. 
In this paper, we finally identify the asymptotic error exponent 
when the optimal strategy for discriminating the states is used. 
A nice feature of such a result is its universality, as it identifies 
the unique metric quantifying the distance of quantum states 
in the i.i.d. setting [5]. 

Distance measures between quantum states have been used 
in a wide variety of applications in quantum information the- 
ory. The most popular such measure seems to be Uhlmann's 
fidelity ((J, which happens to coincide with the quantum 
Chernoff bound when one of the states is pure. The trace dis- 
tance has a more natural operational meaning, but lacks mono- 
tonicity under taking tensor powers of its arguments. The 
problem is that one can easily find states p, a, p', a 1 such that 
Tr|p-a| < Trl/Z-o-'IbutTrl^-a® 2 ! > Ti\p'® 2 - a' m \. 
The quantum Chernoff bound exactly characterizes the expo- 
nent arising in the asymptotic behaviour of the trace distance 
in the case of many identical copies, and therefore does not 
suffer from this problem. Note that a similar situation hap- 
pens in the case of one-copy entanglement versus the asymp- 
totic entanglement entropy. 

In this work we give an upper bound for the probability of 
error for discriminating two arbitrary states. In the particular 
case of a large number of identical copies, this result nicely 
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complements the recent work of Nussbaum and Szkola J3l, 
where a lower bound for the asymptotic error exponent was 
found. These respective upper and lower bounds coincide and 
hence give the exact expression for the error exponent. The 
conjecture of Ogawa and Hayashi concerning the quantum 
Chemoff bound raised in |9] is thus solved. 

Our paper is organized as follows. After the mathematical 
formulation of the problem, we prove a nontrivial and funda- 
mental inequality relating the trace distance to the quantum 
Chemoff bound. Finally, we prove some interesting proper- 
ties of the quantum Chernoff bound and discuss some appli- 
cations. 

PROBLEM FORMULATION 

The optimal error probability of discriminating two quan- 
tum states po and p\ has been identified a long time ago by 
Helstrom |3J]. We consider the two hypotheses Hq and Hi 
that a given quantum system is prepared either in the state po 
or in the state pi, respectively. Since the (quantum) Chernoff 
bound arises in a Bayesian setting, we supply the prior prob- 
abilities 7To and 7Ti, which are positive quantities summing up 
to 1 (the degenerate cases ttq — or — are excluded). 

Physically discriminating between these hypotheses corre- 
sponds to performing a generalised (POVM) measurement on 
the quantum system with two outcomes, and 1 . This POVM 
consists of the two elements {Eq,Ei}, where Eq + E\ = 
1,-Ej > 0. The symmetric distinguishability problem con- 
sists in finding those Eq and E\ that minimise the total error 
probability P e , which is given by 

P e = 7T Tr[Eip ] + 7Ti Tr[E pi}. 

This problem can be solved using some basic linear al- 
gebra. Let us first introduce some basic notations. Abus- 
ing terminology, we will use the term 'positive' for 'pos- 
itive semidefinite' (denoted A > 0) in order to preserve 
trees. We employ the positive semidefinite ordering through- 
out, A > B iff A - B > 0. The absolute value \A\ is defined 
as \A\ := (A* A) 1 / 2 . The Jordan decomposition of a self- 
adjoint operator A is given by A = A + — where A + and 
A_ are the positive and negative part of A, respectively, and 
are defined by. A + := (\A\ +A)/2andA_ := (\A\ - A)/2. 
Both parts are positive by definition, and A + A- = 0. 

Note now that P e can be rewritten as 

P e = 7Ti - Tr[Ei(mpi - ITopo)}- 

This expression has to be minimised over all operators E\ that 
satisfy < E\ < t. The result is that E\ has to be the 
projector on the range of the positive part of (%±pi — ttoPo)- 
We get 

-Pe,min = 7T1 — Tr(7Ti/5i - TT po) + 

= 7ri - (-7T1 - 7r )/2 - Tr \nipi - 7r po|/2 

= 77 (1 - IklPl - T0/>0||l) , 



where 1 1^4| |i = Tr \A\ is the trace norm. 

The basic problem to be solved now is to identify how the 
error probability P e behaves in the asymptotic limit, i.e. when 
one has to discriminate between the hypotheses Hq and Hi 
corresponding to either n copies of po having been produced 
or n copies of p\. To do so, we need to study the quantity 

Pe,n» n ,„:=(lHklP? n -^0P? n |ll)/2- 

It turns out that the behaviour of P e . m i n! „ is exponential 

-P e ,min,n ~ exp {-U^qcb) 

and we will prove that the exponent £qcb is given by the 
following quantity, which can therefore be called the quantum 
Chernoff bound: 

c t 1°§ Pe,min,n 

S.QCB = hm (1) 

n^oo 77 

= - log fmii^ Tr (pV 1 - 8 )^ . (2) 

Note that the quantity Tr (p s cr 1_s ) is well defined and guar- 
anteed to be positive. As should be, this expression for the 
quantum Chernoff bound reduces to the usual definition of 
the classical Chernoff bound £cb when p and a commute: 
for classical distributions po and pi, 




It is truly remarkable that the quantum Chernoff bound is 
given by such a simple expression, looking like an almost 
naive generalisation of the classical Chernoff bound with 
probabilities replaced by noncommuting quantum states. 

The fact that £qc b is lower bounded by the expression on 
the right hand side of (|2ji was proven very recently in 1 8[ (in a 
finite dimensional setting). The fact that this is also an upper 
bound can be inferred from the following theorem, which is 
the main contribution of this paper: 

Theorem 1 Let A and B be positive operators, then for all 
< s < 1, 

Tr^S 1 " 8 ] > Tt[A + B- \A-B\]/2. (4) 

Indeed, let A = nipf n and B = 7r 0/ o^", then the upper 
bound trivially follows from the fact that the logarithm of the 
left hand side of the inequality becomes log(7TQ7ri[~ s ) + 
n log (Tr[pgp} -S ]). Upon dividing by n and taking the limit 
77 — > oo, we obtain the quantum Chernoff bound £qcb, in- 
dependently of the priors tto, tt\ (as long as the priors are not 
degenerate). 

Inequality © is also very interesting from a purely matrix 
analytic point of view, as it relates the trace norm to a mul- 
tiplicative quantity that is highly nontrivial and very useful. 
Note that the optimal measurement to discriminate the two 
sources enforces the use of joint measurements. The particu- 
lar permutational symmetry of A^-copy states however guar- 
antees that the optimal collective measurement can be imple- 
mented efficiently (with a polynomial-size circuit) 1 10[, and 
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hence that the minimum probability of error is achievable with 
reasonable resources. 

PROOF OF THEOREM 1 

Let now move on to prove Theorem^ Note that the proof 
that we present here goes through in infinite dimensions. 

The proof relies on the following technical Lemma, which 
we prove in the Appendix A. 

Lemma 1 Let A, B > 0. Let < t < 1, and let P be the 
projector on the range of (A — B) + . Then 

Tr[PB(A t - B 1 )] > 0. (5) 

Proof of Theorem [7] — We apply Lemma ^ to the case 
t = ,s/(l — s), A = a 1_s and B = b 1 ^ 3 , where a,b are 
positive operators and < s < 1/2. With P the projector on 
the range of (a 1_,s — 6 1_s )+, this yields 

Tr[Pb 1 ' s (a s - b s )} > 0. 

Subtracting both sides from Tr[P(a — 6)] then yields 

Tr^Pia 1 - 3 - < Tr[P(a - b)}. 

Since P is the projector on the range of the positive part of 
(a 1_s - b 1 " 3 ), the LHS can be rewritten as Tr[a s (a 1_s - 
6 1_s ) + ]. Because a s > 0, this is lower bounded by 
Tr^a 1 " 5 - b 1 - 3 )} = Tr[a - a'b 1 ' 3 ]. 

On the other hand, the RHS is upper bounded by Tr[(a — 
6)+]; this is because for any self-adjoint H, Tr[H + ) is the 
maximum of Tr[Qi/] over all self-adjoint projectors Q. We 
thus have 

Tr[a - a'b 1 - 3 } < Tr[(a - b) + ] = Tr[(a - b) + |a - b\}/2. 

Subtracting both sides from Tr[a] finally yields for < 
s < 1/2. The remaining case 1/2 < s < 1 obviously follows 
by interchanging the roles of a and b. □ 

PROPERTIES OF THE QUANTUM CHERNOFF BOUND 

In this Section, we will also study the non-logarithmic vari- 
ety of the quantum Chernoff bound, which we denote here by 
Q{p,(j) := min < s <i Tr^cr 1-5 ]. 

We begin by stating upper and lower bounds on Q in terms 
of the trace norm distance T(p, a) := \\p — ct||i/2; a proof 
can be found in the appendix B. 

1 - Q < T < Vi-Q 2 . (6) 

Based on these bounds, the following properties of the Q- 
quantity and the Chernoff bound can be derived: 

Inverted measure. — The maximum value Q can attain is 
1, and this is reached when p = a. This follows, for example, 



from the upper bound Q 2 + T 2 < 1. The minimal value is 
0, and this is only attained for pairs of orthogonal states, i.e. 
states such that pa = 0. This implies that the Chernoff bound 
is infinite iff the states are orthogonal; this has to be contrasted 
with the asymmetric error exponents occuring in the context 
of relative entropy, where infinite values are obtained when- 
ever the states have a different support. 

Convexity in s. — The function to be minimised in Q is 
s i— > Ti[p s cr 1_s ]. It is important to realise that this function is 
convex in s £ [0, 1], because that means that the minimisation 
has only one local minimum and therefore this local minimum 
is automatically the global minimum. This is an important 
benefit in actual calculations. 

Indeed, the function s t— * a; s y 1_;i is convex for positive 
scalars x and y, as one easily confirms by calculating the sec- 
ond derivative x s y 1 ~ s (logx — logy) 2 , which is non-negative. 
Consider then a basis in which p is diagonal and given by 
p = Diag(Ai, A2, . . .). Let the eigenvalue decomposition of a 
(in that basis) be given by a = U Diag(/ii, /12, . . .)£/*, where 
U is a unitary. Then Tr^o- 1 - 5 ] = Kp)~ s \U l3 \ 2 . As 
this is a sum with positive weights of convex terms X s p^~ s , 
the sum itself is also convex. 

Joint concavity in (p, a). — By Lieb's theorem flUl . 
Tr[p s a 1 ^ 11 ] is jointly concave in (p, a). Since the quan- 
tum Chernoff bound is the pointwise minimum of Tr[/9 S CT 1_S ] 
(over a fixed set, namely over s S [0, 1]), it is itself jointly 
concave as well. The Chernoff bound is therefore jointly con- 
vex, just like the relative entropy. 

Monotonicity under CPT maps. — From the joint concavity 
one easily derives the following monotonicity property: for 
any completely positive trace preserving (CPT) map $, 

Q($(p),$( ( 7))>Q(p,a). (7) 

To prove this, one first notes that Q is basis independent, i.e. 
is invariant under unitary conjugations: 

Q{UpU*,UaU*)=Q(p,a). 

Secondly, Q is invariant under addition of an ancilla: let r be 
the (normalised) ancilla state, then 

Q{p ® t, a <gi t) = Q(p, cr); 

this is because Tr[(p ® r) s (cr ® r) 1 " 8 ] = Tr^cr 1 ^] Tr[r]. 
Exploiting the Stinespring form of a CPT map, the monotonic- 
ity statement follows for general CPT maps if we can prove it 
for the partial trace map. As noted by Uhlmann f^fol . the 
partial trace map can be written as a convex combination of 
certain unitary conjugations. Monotonicity of Q under the 
partial trace then follows directly from its concavity and its 
unitary invariance. 

Continuity. — By the lower bound Q + T > 1, 1 — Q is 
continuous in the sense that states that are close in trace norm 
distance are also close in 1 — Q distance: < 1 — Q < T. 

Relation to Fidelity — If one of the states is pure, then Q 
equals the Uhlmann fidelity. Indeed, assume that p\ = \ip)(ip\ 
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is pure, then the minimum of the expression Tr(p|/9 2 ~ s ) i s 
obtained for s = and reduces to (ip\p2\ip). As shown in the 
appendix B, the fidelity is always an upper bound to Q. 

Relation to the relative entropy — Just as in the classical 
case, there is a nice connection between the quantum relative 
entropy and the Chernoff bound. By differentiating the ex- 
pression Tr(p s cr 1_s ) with relation to s, one observes that the 
minimum (which is unique due to convexity) is obtained when 

Tr(pV~ s logp) = TrGoV- 8 logcr). 

By the cyclicity of the trace, one easily verifies that this is 
equivalent to the condition that 

S(r s \\p) = S(r s \\a) 

with the quantum relative entropy Tr(Alog^4 — 

A log B) and r s defined as 

s Tr J o s cr 1_s 

Note that t s is not a state, because it is not even self-adjoint 
(except in the commuting case). Nevertheless, as it is basically 
the product of two positive operators, it has positive spectrum, 
and its entropy and the relative entropies used in (jSJi are well- 
defined. The value of s for which both relative entropies co- 
incide is the optimal value s* . This t s * can be considered the 
quantum generalisation of the Hellinger arc and interpolates 
between two different quantum states, albeit in a rather spe- 
cial (unphysical) way. 

Metric. The quantum Chernoff bound (or its non- 
logarithmic variety) between two infinitesimally close states 
p and p — dp induces a metric that gives a geometrical struc- 
ture to the state space. In Appendix C this metric is shown to 
be 



ds 2 



1- min Tr[p s (p- dp) 1 

0<s<l 



i v l(W)l 2 



where p = J^. Aj|i)(i| is the eigenvalue decomposition of p. 



CONCLUSION 

We have identified the exact expression of the quantum 
generalization of the Chernoff bound, which allows to quan- 
tify the asymptotic behaviour of the error in the context 
of Bayesian discrimination of different sources of quantum 
states. This resolves a long-standing open question. Our main 
theorem (Theorem 1), which gives a computable lower bound 
to the trace norm difference of two states in the many-copy 
regime, may also find other relevant applications in and out- 
side the field of state discrimination. 
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APPENDIX A 

Proof of LemmaUi — We exploit the integral representation 



t sin(i7r) / ' n.r 



dx 



a + x 



(9) 



which is valid for a > and < t < 1 1 14]. Extending this 
representation to positive operators in the usual way, we get 

Tr[PB(A< - B<)] = [ +OC dxx 1 - 1 



Tr[PB(A(A + x)- 1 - B(B + x)- 1 )] 

If we can prove that the integrand is positive for all x > 0, then 
the integral itself is also positive. To do so, we first reduce 
the integrand to yet another integral and then prove that the 
integrand of that integral is positive. Let A = A — B. Now 
note 

A(A + X)- 1 - BiB + x)- 1 
= (B + A){B + A + x)- 1 - B(B + X)- 1 

r 1 d 

= / dt — (B + tA)(B + tA + x)- 1 
Jo dt 

= / dt x{B + tA + x)- 1 A(B + tA + x)' 1 . 
Jo 

Here, the last equality can be shown as follows: with B' := 
B + tA, 



-(B + tA^B + tA + x)- 1 
dt 

= A{B' + X)- 1 -B'(B' + x)- 1 A{B l + x)~ x 



x(B' + x y L A(B' +x) 



Therefore, 
Tr[PB(A(A 



-BiB + x)- 1 )} 



= x [ dt Tr[PB(B + tA + x)- 1 A(B + tA + x)- 1 }, 
Jo 

where we note that B + t A is positive for < t < 1. 
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If we can show that the integrand is positive for < t < 1, 
then the integral itself is also positive. In the following we 
absorb t in A and write C as a shorthand for B + t A. What we 
have to prove then is that for all B, C > 0, a; a non-negative 
scalar, and P the projector on the range of (C — B) + , 

Tr[PB(C + x)- 1 A(C + X)- 1 } > 0. (10) 

Let A have the Jordan decomposition A = A + — A_ . Thus P 
is the projector on the range of A + . We introduce the symbol 
V = (C + x)^ 1 > 0. We choose a basis in which A and P 
can be partitioned as 

*-(V-L> *-(!2> 

In that same basis, V can be partitioned as 

v= (V n V 12 \ 
\ V 2 i V22 J ' 

For ease of notation, the subscript ij will henceforth refer to 
the (i,j)-th block of an operator valued expression. Inequality 
J lOt can then be rewritten as 

Tr[P(C - A)(VAV)] = Tr[(C - A)(VAV)]n > 0. 

Noting that C — A = V^ 1 - A — x, the LHS is equal to 

TrKF" 1 - A - x)(VAV)]u 
= Tr[AV- {A + x)(VAV)] n 
= Tr[A+Vu - (A+ + x)(VAV)u] 
= Tx\A+{V -VAV^-xiVAV)^. (11) 

Because of the positivity of B, we have V^ 1 > A + x, which 
implies V = VV- l V > V{A + x)V = VAV + xV 2 . As the 
diagonal blocks of a positive operator are themselves positive, 
this further implies 

V u - (VAV)u > x(V 2 ) n . 

Inserting this in Q 1> gives 

Tr^y- 1 - A-a;)(yAF)] u 
= Tr[A+(V- VAV)u -x(VAV) n ] 
> TrlA+a;^ 2 )!! - x(VAV) n ] 
= xTriA+iV^n-iVAV)^ 

= xTi[A + {V n Vu + V 12 V 2 i) - {VnA+Vu - V 12 A_V 21 )] 
= xTt[A + Vi 2 V 2 i + V 12 A^V 21 }. 

By the fact that V12 and V 2 \ are each other's adjoint, the latter 
expression is positive, which finally proves the statement of 
the Lemma. □ 

APPENDIX B 

Theorem^ applied to normalised states immediately gives 
the lower bound 

Q + T>1. (12) 



Below we provide an upper bound on Q that is valid for any 
pair of states. 

By definition, for any fixed value of s between and 1, the 
quantity Tr [p s (T 1_s ] is an upper bound on Q. In what follows, 
we set s — 1/2. Furthermore, by replacing the trace with the 
trace norm, we get an even higher upper bound. Indeed, 

Q < Trip 1 / 2 * 1 ' 2 } 

= i^v/y^iii 

< \\p 1/2 v 1/2 \\i- (B) 

In the last line we have used the fact (|14], Prop. IX. 1.1) 
that for any unitarily invariant norm |||AB||| < |||BA||| if 
AB is normal. In particular, consider the trace norm, with 
A = p 1 / 4 * 1 / 2 and B = p 1 / 4 ; then AB is self-adjoint, hence 
normal. 

What we obtain as the RHS is the so-called Uhlmann fi- 
delity F between the states p and a: 

F := Tr[{p 1/2 ap 1 / 2 ) 1 / 2 ] = Wp 1 ' 2 * 1 ' 2 ^. 

This quantity can be regarded as the generalisation of the 
Bhattacharyya coefficient to the quantum case. For classical 
distributions po and pi, the Bhattacharyya coefficient is de- 
fined as |1J B(p ,pi) := £V \/po(«)pi(«)- We have thus 
just shown that F is an upper bound to Q. 

Furthermore, by a result of Fuchs and van de Graaf 1 15], 
1 — ^F is lower bounded by the square of the trace distance 
T: 

T 2 < 1 - VF. (14) 
Combining this with inequality dl 31 yields the upper bound 

Q 2 +T 2 <1. (15) 

There is a nice direct proof of the latter inequality that cir- 
cumvents the proof of J14i and goes through in infinite di- 
mensions. We state it in terms of general positive operators: 

Theorem 2 For positive operators A and B, 

||A-5||2+4(Tr[A 1 / 2 B 1 / 2 ]) 2 < (Tr(A + B)) 2 . (16) 

Proof. Consider two general operators P and Q, and define 
their sum and difference as S = P + Q and D = P — Q. We 
thus have P = (S + D)/2 and Q = (S - D) /2. Consider the 
quantity 

PP*-QQ* - ^((S + D){S + D)* - (5 -£>)(5 -£>)*) 

= \{SD* +DS*). 
Its trace norm is upper bounded as 

||S.D*+DS*||i/2 < (\\SD*\\ l + \\DS*\\ l )/2 
= \\SD*\\ X 
< \\S\\2\\D\\ 2 . 
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In the last line we have used a specific instance of Holder's 
inequality for the trace norm ([14] Cor. IV.2.6). Now put 
P = A 1 / 2 and Q = B 1 / 2 , which exist by positivity of A 
and B, and which are by themselves positive operators. We 

get S,D = A 1 / 2 ± B 1 / 2 , hence 

P-B||i < WA 1 / 2 + B^ 2 \\ 2 HA 1 / 2 --B 1/2 || 2 , 

which upon squaring becomes 

\\A-B\\\ < Tt{A 1 ' 2 + B 1 ' 2 ) 2 Tr{A 1 ' 2 - B 1 ' 2 ) 2 
= Ti(A + B + A 1/2 B 1/2 + B 1/2 A 1/2 ) 

x Tr(A + B- A 1 I 2 B 1 ' 2 - B 1 i 2 A 1 ' 2 ) 
= (Tr(A + B) + 2Tr(A 1/2 B 1/2 )) 

x (Tr(A + B)-2 Tr(A^ 2 B^ 2 )) 
= (Tr(A + B)) 2 -4(Tr(A 1/2 B 1/2 )) 2 . 

□ 



Together with the lower bound Q + T > 1 we can now 
bracket the trace distance in function of the Q quantity: 



1 - Q < T < ^l-Q 2 < 1 - Q 2 /2, 



(17) 



where the last inequality becomes a very good approximation 
for small values of Q. 



APPENDIX C 

Derivation of the quantum Chernoff metric: The goal is to 
calculate 



ds = 1 — min Tr\p s (p — dp) 

Q<s<l 



(18) 



Here we will use the integral representation l|9) as well as 
its derivative, 



, , sin(tn) f + °° , x 
to,*' 1 = — — / dx 



{a + xf 



(19) 



which holds for a > and — 1 < t < 1. 

In particular, using (|9} and the convergent sequence 

l/(a - b) = a^ 1 + a~ 1 ba~ 1 + a~ 1 ba~ 1 ba~ 1 + . 

one can write (p — dp) 1 ^ 8 up to second order in dp, 



X 



{p-dpf s = c s I dx (p-dp) 

J a p- dp + x 

f + °° _ ( 1 

~ c s / dx x s (p — dp) 

Jo \P + x 

1 j 1 1,1,1 
H dp 1 dp dp 

p + x p + x p + x p + x p + x 

where c s = ir^ 1 sin[s7r]. 



Inserting this expansion in Jl 8I > one finds 



ds 2 



0<s<l 



dx Tr 



{p + x)' 



; P S dp 



(p + x} 



; P S dp 



p + x 



dp 



(20) 



The first term in the integrand vanishes, as can be seen by 
using d!9t and Tr dp = 0, while the second term can be com- 
puted in the basis that diagonalizes p = J^i 



ds 2 = max y c s / dxx 

0<s<l^ ./ n 



A? 



(Aj +x) 2 {\ 3 +x) 



\(i\dp\j)f 



— max 



\(i\dp\j)\' 



c - i ^-y (Ai - X 3 ) 2 



(Ai + Xj - AfAj s - X s j X l s ) 



v J 
l.v- \(i\d P \j)\ 2 



2 V^ + V A i) 



2 ' 



(21) 



where in the second equality we have used dp — dp^ and 
symmetrized the factor that multiplies | (i\dp\j) \ 2 in the sum. 
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